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Introduction 

Recent research has revealed a surge in the application of Stochastic Differential 
Equations (SDEs) in the modeling of infectious diseases. Factors emanating from this surge 
has been the ability of stochastic differential equations to lay one of the most imperative 
fundamental differences in the asymptotic dynamics of the deterministic and stochastic 
epidemic models [1]. These fundamental differences include the convergence of the sample 
path of stochastic models to the disease-free state and that of the deterministic model to 
the endemic equilibrium. Stochastic phenomenon occurs naturally in our environment, and 
when small number of reacting molecules is involved in the modeling system, deterministic 
models become inappropriate [2]. Stochastic models have properties which include the 
probability of an outbreak, the quasistationary probability distribution, the final size 
distribution of an epidemic and the expected duration of an epidemic that makes it unique 
from others [3]. Stochastic Processes can have moments and covariance function associated 
with them, which are functions of time, as do random variables [4]. 

Stochastic epidemic models have been studied extensively, and the outcomes have 
been astonishing as they have provided an alternative means of investigating epidemic 
outbreaks in a small populations coupled with better elucidating of the transmission 
dynamics of epidemic outbreaks in the early stages of an epidemic outbreak, insinuating the 
robustness of the deterministic model with respect to stochastic perturbations [5]. Kiouach 
and Omari [6], considered an epidemic model of type SIRS and the stability of the interior 
equilibrium point was discussed. The model was modified to stochastic version and the 
stability properties of the stochastic model were investigated. Maroufy et al [7], studied a 
classical model of a SIRS epidemic in an open population. They deduced an explicit formula 
by which the lower bound of the density of the infective can be computed. The Model was 
then extended to stochastic version and Lyapunov functional was employed to investigate 
the global stability of both the deterministic and the stochastic model. Lahrouz et al [8], 
presented an epidemic SIRS model with saturated incidence rate and disease-inflicted 
mortality. The model was modified into stochastic version and the global stability of the 
models were studied using the tool of Lyapunov function. Further, the global existence and 
positivity of the solution was verified as well as the global stability in probability proved 
under suitable condition by the perturbation of white niose. Murkherjee [9], considered 
deterministic and stochastic models for prey-predator system where the prey population 
is infected by micro parasite. The stochastic stability properties of the models were 
investigated which suggested the robustness of the deterministic model to stochastic 
fluctuations. Clancy [10], formulated a stochastic SIS epidemic model for transmission of 
infectious disease through a population, considering direct host-host transmission and 
indirect transmission through free-living infectious stages. Guoting and Tiecheng [11], 
investigated a stochastic version of the SIR model. The stability in probability of the steady 
state of the system was proved under favorable conditions with white noise perturbations. 
Mukherjee et al [12], studied the behaviour of a plant herbivore for both discrete and 
continuous model with stochastic perturbation. The existence and stability of two fixed 
points were investigated which suggested local asymptotic stability in probability for the 
stochastic model for certain strengths of white noise. Bifurcation diagrams and time series 
plots were obtained for the model. 

Entrenched in these works are devotees of outstanding achievements and contributions 
towards scientific community on the mathematical modeling of infectious diseases. In this 
section, the dynamics of deterministic model of malaria disease and its stochastic variant 
are presented. The stochastic version of the model is derived by introducing random 
perturbation to the equations of system (1). Numerical examples are provided to illustrate 
the explicit differences between the deterministic model and the stochastic system. 

The Stochastic Epidemic Model Formulation 

The stochastic model of the deterministic system (1) is obtained by introducing 
stochastic perturbations in the main parameters of the deterministic system. Thus we 
permit stochastic perturbations of the variables S H , E H , I H , R H , S y , E y , I y around 
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their values at positive equilibrium E* [13,14,15]. Further, we assume 
that the white noise of the stochastic perturbations of the variable 
around values of E* are proportional to the distances s h ,e h ,i h ,r h ,s v ,e v j v 
from s* h ,e* h ,i*„,RhXXX' Hence, the stochastic version of model (1) is 
given by 


dS it — 


n H + Ph R H ^vPh j.j R H Ph R H 


dt + Q { (S H -S H )dB 1 


dE H 


°vPh —Sh 

v H N„ H 


(v h + p H )E h 


dt + Q 2 (E h -E* h )dB 2 


dI H ~ \y h^h (Yh Ph )^h ] dt + Q 3 (I h I H )dB 3 
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dl v = [v,E y -/j. v I v \dt + Q 7 ( I v -ry)dB 7 


dt + Q 5 (S v - S* v )dB 5 
dt + Q 6 (E v - E* v )dB 6 


With Q., for = 1, 2.7 are real constants and B = 1, 2, 3, 4, 5, 

6, 7 are independent wiener processes. 

Numerical Method for SDEs-Euler-Maruyama 

This section focuses on stochastic modeling application and 
numerical simulations of trajectories. The work employs the 
fundamental works which allow the development of effective 
numerical methods for stochastic differential equations (SDEs), as 
well as providing the basic information necessary to understand the 
nature of SDEs and how to stimulate them by using the tool of Euler- 
Maruyama method [16,17,18]. The method has been proved to have 
a strong global convergence under the assumption that the stochastic 
differential equation is locally Lipschitz and that the moment of the 
exact and numerical solution are bounded for some p > 2. 

Consider a scaler, autonomous SDE in integral form 

X(t) = X 0 + £ f(s,x(s))ds + £ g(s,x(s))dw(s), 0 < t < T (2) 

Where f and g are scaler functions and the initial condition x Q is a 
random variable. 

If x(t) is the solution to (2), then the solution x(t) is a random 
variable that arises when we take the zero step size limit in the 
numerical sense. Hence, the differential equation form of (2) can be 
written as 

dx(t) = f (/, x{t))dt + g(t, x(t))dw{t ), with X(0)=0 < t < T (3) 

Now, from equation (3), if g= 0 and x Q is constant, then the 
problem becomes deterministic and (3) reduces to autonomous 
ordinary differential equation 

~~ = fit, x(t)), with X(0)= x 0 (4) 

Applying Euler-Maruyama method to (3) over (0, T), we first 

discretize the interval by letting At = j 

for some positive integer L and T.= j At . Hence, our numerical 
approximation to X (T.) will be denoted by X,. Therefore, the Euler- 
Maruyama (EM) method becomes 

Xj =X j _ l +f(X j yAt + g(X j _ l )(W(T j )-W(T j _ 1 )),forrl, 2 , .L ( 5 ) 

Equation (5) arises from the integral form 

X(Tj ) = X(Tj_,) + ££ f(s, x(s))ds +££ g(s, x(s))dw(s) ( 6 ) 

Further, equation (6) reduces to Euler's method when g = 0 and 
x 0 is constant (see [19,20,21]). 

Numerical Examples and Discussion 

Here, we illustrate with figures the dynamics of the deterministic 
and stochastic systems, and gives an explicit difference in the models 


by carrying out numerical simulation for hypothetical set of parameter 
values. To demonstrate the differences, we simulate the stochastic 
and the deterministic systems by using the set of parameter values in 
table 1. To confirm the deterministic plots, we choose arbitrary white 
noises Q 2 , Q 3 , Q 4 , Q 5 , Q 6 and Q 7 of equal strength 0.5, and shows the 
fluctuations in the trajectories of the stochastic system. We can see 
that the trajectories of the stochastic plots displayed are the same as 
the trajectories of the deterministic model during a finite time frame. 
This suggest that even though fluctuation occurs, they result in the 
same behavior as the deterministic one. The deterministic approach 
has limitations that the stochastic handles in a more realistic way. 
Thus, the deterministic model gives the same results any time the 
simulation is run with the same parameter and initial values. This 
might be mathematically correct, however, this is not the case in real 
epidemic situation, as there may be many parameters which we cannot 
not model entirely realistically. Modelling them deterministically, we 
lose some of the complexity of the system. Hence, it is appropriate 
to assume a stochastic perturbation. Hence, the differences in the 
dynamics of the models are given by the figures 1-12. 

Conclusion 

In this section, the dynamics of the deterministic malaria model 
with standard incidence rate and its stochastic version are presented. 
The stochastic version of the model was constructed by employing 
stochastic fluctuations in the main parameter of the deterministic 
model. Our main purpose of the study was to investigate the dynamics 
of the deterministic and stochastic systems, by giving numerical 
examples to illustrate the explicit difference in the models. The 
numerical simulations for the models show that the traj ectories of the 
stochastic plots were the same as the trajectories of the deterministic 
model. Further, from the stochastic plots, the simulations show 


Table 1: Description of variables and parameters of the Malaria 
Model (1) 


Parameter 

Description 

Estimated 

value 

Reference 

n H 

Recruitment rate of human 

2.5 

[9] 

Pn 

Transmission probability from 
an infectious 

mosquito to a susceptible 
human 

0.9 

[9] 

A. 

Transmission probability from 
an infectious 
human to a susceptible 
mosquito 

0.8 

[9] 

PvH 

Transmission probability from a 
recovered 

human to a susceptible 
mosquito 

0.009 

[9] 

Ph 

Natural death rate of human 

0.00004 

[9] 

S H 

Death rate from the disease 

0.00354 

[9] 

Yh 

Infectious human recovery rate 

0.003704 

[9] 

v H 

Progression rate from EH to IV 
class 

0.08333 

[9] 

V v 

Progression rate from EV to IV 
class 

0.1 

[9] 

Ph 

Rate of loss immunity in human 

0.0146 

[9] 

<7 V 

Recruitment rate of mosquitoes 

500 

[9] 

(J v 

Biting rate of vector 

2.9 

[9] 

Pv 

Natural death rate of vector 

0.0714 

[9] 

a H ,a v 

The level at which the force of 
infection saturates 

(0,1) 

[9] 
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Figure 3: The plot represents population of Exposed Human 
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Figure 8: Stochastic plot of Susceptible Mosquito 
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Figure 10: Stochastic plot of Exposed Mosquito 
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Figure 12: Stochastic plot of Infected Mosquito 


an initial random fluctuation of the stochastic trajectories, until 
they eventually approach asymptotic level. The numerical example 
shows that the stochastic system is globally asymptotically stable in 
probability when the densities of white noise are less than certain 
threshold parameters. However, if these densities of white noise are 
zero, it means there are no stochastic environmental factors on the 
population and hence no stochastic perturbation. 
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